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^ Abstract. We analyze optical spectra of the high temperature superconductor HgBa2Cu04+5 

^ using a minimal model of electrons coupled to bosons. We consider the marginal Fermi liquid 

^ ^ theory and the spin fluctuation theory, as well as a histogram representation of the bosonic 

spectral density. We find that the two theories can both be used to describe the experimental 
!/3 data provided that we allow for an additional scattering channel with an energy of 55 meV. 

B 

^ 1. Introduction 

O The puzzle of superconductivity in the cuprates remains unsolved. Over the years it has become 

, clear that the nature of the normal state is an equally important puzzle to solve. It now seems 

established that the strongly overdoped cuprates posses a well defined Fermi surface [H ?,|2l ?] 
indicating that at least in this part of the phase diagram at low energy a Fermi liquid picture 
might apply for the normal state. The underdoped cuprates are characterized by so-called Fermi 
arcs: incomplete pieces of Fermi surface [?]. One might expect that in this part of the phase 
diagram the Fermi liquid paradigm is no longer appropriate and the physics is indeed described 
by that of a lightly doped anti- ferromagnetic insulator. Recent de Haas - van Alphen experiments 
however seem to indicate that the arcs should be interpreted in terms of small pockets [3] thus 
placing the state underlying the pseudogap state also in the Fermi liquid regime. Regardless the 
nature of the superconducting and normal state, the transition from one to the other is driven 
sjjl by the minimization of the free energy, which consists of potential {Epot) and kinetic energy 

•'^ (Ekin). In conventional superconductors the transition is driven by a lowering of the potential 

rN energy while at the same time the formation of the Cooper pairs goes at the cost of a smaller 

amount of kinetic energy. Several experiments [H El 1^ have found that this situation might be 
different in the cuprates. In |3] the experimental results could be explained if it was assumed 
that superconductivity is driven by a gain in the c-axis E^in while in the latter two cases the 
result can be interpreted in terms of a gain in a6-plane E^in- In both cases it was argued that the 
different sign of E^in with respect to conventional superconductors arises from the properties 
of the normal state (see i.e. Ref. [?] for the first case and Ref. [?, ?] for the second). In 
later experiments O El |9] the dependence of the change in E^in on doping and the number of 
layers was studied. The results of these studies are summarized in figure [T} The quantity AW 
is a measure of the kinetic energy change at the superconducting transition as defined in |10] . 
Interestingly, AW changes sign for slightly overdoped samples. Assuming that the mechanism 
of superconductivity does not depend on doping, the sign change must arise from a change in 
the properties of the normal state. 
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Figure 1. Doping dependence of the 
spectral weiglit cliange AW at Tc- The 
numbers are taken from Ref. O [7] for Bi- 
2212, from Ref. [8] for Bi-2223 and from 
Ref. [9 for Hg-1201. The result for Nd-123 is 
unpublished. This result is similar to the one 
obtained by the group of Bontemps [?]. 



A possible interpretation of this result is that the superconducting dome covers a quantum 
critical point [?, |Tl]. In this scenario the underdoped and overdoped materials would be 
distinct phases explaining the change in Fermi surface topology and the change in sign of 
the superconductivity induced spectral weight. It has also been suggested that the optical 
properties of optimally doped Bi-2212 [TT and Bi-2223 [?] can be understood from quantum 
critical fluctuations. In [12^ it was pointed out that at low energy the optical spectra could 
be described as a universal function of w/T while in the mid-infrared range the conductivity 
could be described by a powerlaw instead. However, as pointed out in [?] these two observations 
cannot be made consistent with each other without assuming a non-universal background to the 
optical conductivity. 

Norman and Chubukov [13] showed that the powerlaw behavior could also have another 
interpretation: it follows that the optical conductivity has an approximate powerlaw behavior if 
one assumes a model in which electrons interact with a broad bosonic spectrum (n(u;)). Using 
the same model the anomalously large temperature dependence of the normal state spectral 
weight could also be explained [H] but not the observed change in kinetic energy at Tc [?]. 
Stimulated by these positive results we analyzed the optical spectra of 10 different compounds 
using a histogram representation for the bosonic spectrum [?] and found that this spectrum 
strongly depends on doping and temperature. The doping dependence of n(u;) might give an 
explanation for the very different Fermi surfaces observed by ARPES for under- and overdoped 
cuprates, while the temperature dependence of the spectra could give an explanation for the 
(approximate) uj/T scaling observed in [?]. Here we will discuss two possible models for the 
normal state and compare it to the earlier obtained histogram representation. 

2. Methods 

In normal metals the electron-phonon (EP) interaction can be described in a framework 
where fermionic quasiparticles, the electrons, interact with a spectrum of bosonic modes. 



the phonons. The function describing the energy dependent couphng is often indicated as 
a'^F{u)). A closely related theory is the Eliashberg theory of superconductivity describing the 
occurrence of superconductivity in normal metals. It is an improvement of the BCS theory of 
superconductivity because it takes into account the retardation of the EP interaction. a'^F{uj) 
also appears in this theory and is the main determining factor of the critical temperature at 
which superconductivity occurs. It is also possible to consider bosons other than phonons, for 
example plasmons, magnons or excitons. An important distinction between these bosons and 
phonons is that they are made up out of the same degrees of freedom which form the Cooper 
pairs: the electrons themselves. 

We compare different models for the electron-boson coupling function, indicated here by the 
symbol n(u;) . The calculation of the optical conductivity for given n(u;) is straightforward (see 



Appendix A). The inverse problem of extracting H{lu) from experimental data can be done, but 



with limited accuracy due to the progression of errors due to experimental noise [151 IIS] • Several 
theoretical models for the cuprates predict that the electrons interact with a bosonic spectrum. 



2.1. Marginal Fermi liquid. 

The Marginal Fermi Liquid (MFL) theory |17| [TT] predicts that n(ci;) is given by, 

nA/FiM = Atanh(^)/(w,Wc,A) (1) 

where A is an overall coupling constant, T is the temperature and /(u;,a;c,A) is a high 
energy cutoff function the precise shape of which is unimportant. We have used /(u;,Wc, A) = 
1/(1 + exp{ ^^^""•^ )) as proposed originally by Varma [TT]- An important feature of this spectrum 



is that it only depends onuj/T. It is easy to see using Eq.'s A.l A. 4 that in this case one retrieves 
the linear temperature dependence of the resistivity seen in many experiments. It would also 
give rise to a conductivity that scales as a{uj,T) = T~^g{uj/T) as proposed in [?]. 



2.2. Spin fluctuation theory. 

Several other theoretical models are based around the notion that anti-ferromagnetic fluctuations 
play an important role in understanding the physics of the (underdoped) cuprates [?, 18, ?]. In 
these models one uses the imaginary part of the dynamic spin susceptibility y^'{oj,q = 0) as a 
measure for the bosonic spectral density which in the ungapped state is given by, 

^i-SFV^) = n W < (2) 

7^ -|- oj'^ 

Here T is a coupling strength and lOc a cutoff. The spectrum has a maximum determined by 
a; = 7. We refer to this model as the spin fluctuation (SF) model. Note that this form is 
valid only if there is no (pseudo)gap. Although the temperature dependence is not explicitly 
mentioned in Eq. [2] the parameters 7 and V depend on temperature. lispiyj) does not scale as 
u}/T unless both 7 and T are depending linear on temperature. The temperature dependence 
of the MFL and SF models are an important point of contrast with the EP interaction which 
leads to a temperature independent a^F{uj). The fact that most experimental studies find that 
for the cuprates n(u;) has to be made temperature dependent is one of the arguments against 
the interpretation of n(a;) in terms of phonons. 



2.3. MFL and SF + phonon models 

Finally we allow for a Lorentzian oscillator in addition to the MFL and SF model given by. 



Here /i is the oscillator strength, Ti is the width and uji is the center frequency. Such a peak 
can be used to describe the coupling of electrons to phonons or to the spin resonance. 



2.4- Histogram representation. 

Although we cannot directly invert the experimental data to obtain n(a;) we can use a Levenberg- 
Marquardt optimization routine to determine the important features in n(w) . We do this by 
making a histogram of Il{uj) using i blocks with flexible widths and heights with a total of 2i 
parameters, 

^hg{^) = fi < w < cjj (4) 

where lvq = and fi is the height of the ith block. We find that i=6 is the optimal number 
of blocks giving sufficient detail to n(a;) without overdetermining it. For ujq < to < coi we use 



^hg{^) = fi^ to circumvent the divergence of the integral in Eq. A. 2 The resulting histogram 
allows us to extract a rough impression of n(u;) and give an indication of the important features 
in the spectrum. This model function allows us to optimize n(a;) such that adding more detail 
does not significantly improve the result. This model was used in a recent study of 10 different 
compounds [?]. We showed that II}{g{^) has three important features: 

• tlHci^) has to be temperature dependent. 

• IIhg{^) extends to very high energies (about 400 meV) as compared to standard electron- 
phonon systems 

• HHci^) is dominated by a peak (55 meV) in the phonon range. 

The novel idea that follows from these observations is that the glue function of the cuprates is 
made up out of electronic and phononic contributions. Below we will come back to this point 
and try to estimate the separate contributions. 

3. Results 

We analyze optical spectra of optimally doped HgBa2Cu04+5 (Hg-1201). This high temperature 
superconductor presents us with a good opportunity to test the above models without too many 
complications. It is a single layer compound with a simple tetragonal structure and it has a 
high critical temperature of Tc ~ 97 K. The optical data has been presented before in Ref. 
[S]. We do not take into account the opening of a (pseudo-)gap in our analysis so we can only 
apply it to temperatures of 100 K and higher. Recent neutron scattering experiments [?] show 
that at optimal doping the pseudogap temperature is probably smaller or equal to the critical 
temperature for these materials at optimal doping, although specific heat measurements on the 
particular sample used in this study show an onset of the specific heat jump at slightly higher 
temperature [9]. In figure [2] we show the reflectivity for Hg-1201 together with the optimized 
model curves obtained with Eq.'s ([T|, ^ and (Q. The corresponding spectral functions are 
shown in figure |3j Figure |2^ and show the results for the MFL model. At 290 K we 
can make a reasonable fit to the reflectivity but with decreasing temperature the fit becomes 
progressively worse. Note that we try to fit the overall frequency dependence of R{uj) and not 
the sharper phonon structures below 80 meV. It is important to note that in order to optimize 
the mean squared deviation, (see appendix), we have to make the coupling constant A in Eq. 
[T] temperature dependent. This means that already at this level the uj/T scaling is lost even 
for the MFL model. Our result obtained at room temperature is similar to that obtained in 
an earlier study [?] . We find that if we add an extra scattering channel to the MFL model the 
fit is greatly improved (see below). I[sf{^) (panels [ija and jsj)) shows a similar trend as seen 
for IIm^{uj) : a reasonable fit at high temperature but less good at low temperature. Finally, 
andjsj: show I[hg{^) which gives the best description of the optical data. The panels 
ilgive a lot of information on the structure of Il{co) . As mentioned above it shows three 
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Figure 2. Reflectivity of Hg-1201 together Figure 3. Temperature dependent spectral 

with fits for three different models, (a): functions corresponding to the fits shown in 

MFL model, (b): SF model and (c): the figure [2} (a): MFL model, (b): SF model 

histogram. and (c): histogram. 



important features: (i) n(a;) has to be made temperature dependent, (ii) there is a high energy 
scale that is determined by the cutoff in 11(0;) around 300 meV. (iii) there is a low energy scale 
determined by a maximum in n(u;) around 55 meV. We do not find evidence for ujjT scaling of 
the function n(Lij) from the above three model calculations. Even in the case of the MFL model 
we have to adjust the coupling constant A for each temperature and this destroys the perfect 

w/T scaling. 

To objectively compare the different models we use x^i (Eq. A. 5). Table[l]lists the optimized 
for the fits shown in figure [2} From the table it follows that overall the SF model performs 
better than the MFL model. The reason for this is that the SF model has a maximum around 55 
meV similar to the peak seen in the histogram representation. The MFL model in contrast has 
a broad maximum around 150 meV. The histogram representation suggests that n(u;) consists 
of two contributions: a peak on top of a broader continuum. In [?] we showed that this is true 
for a whole range of dopings and temperatures of single, double and triple layer cuprates. It is 
therefore tempting to ascribe the peak to oxygen vibrations present in this energy range while 
the continuum indicates coupling to either spin fluctuations or loop current excitations. There 
are different interpretations of our observations possible however. 

A first possibility is that one of the assumptions underlying the strong coupling theory is 
simply no longer valid. In these materials the Fermi energy is quite small while the average 
mode energy is relatively high, which could mean that the Migdal approximation, which consists 
of neglecting vertex corrections, is no longer valid. This could mean that the high energy part of 
the background would be absent if these corrections are properly taken into account. However, 



Table 1. defined by Eq. A. 5 for the three different models. The values in parentheses are 
obtained when a Lorentzian peak (Eq. pi) is added to the spectral function. 



Model 290 K 



200 K 



100 K 



MFL 117 (13) 127 (21) 163 (22) 
SF 23 (15) 37 (21) 74 (18) 
HG 11 19 13 



Chubukov and Schmalian [20] considered 3D fermions coupled to a massless bosons spectrum 
and showed that vertex corrections can be neglected even in the strong coupling limit. 

A second possibility is that we are interpreting parts of the electronic response in terms of 
n(a;) which is actually due to strong correlation effects. According to Anderson [21] the physics 
is determined by two energy scales: the antiferromagnetic exchange coupling J and the Hubbard 
repulsion U. Since J ~ 0.15 eV it is possible that we are incorporating features in Pi{uj) which 
are actually related to RVB like fluctuations related to J. 

The comparison of the panels in figure [3] suggests that we may get a better fit if we add a 
narrow peak to the MFL or SF model. As shown in table[T]the addition of this peak dramatically 
improves the quality of the fit, in particular for the MFL model. We can now easily separate 
the contribution due to the peak and the one due to the continuum which is not possible for 
the histogram representation. For each spectrum we can calculate the coupling constant which 
is given by, 



and separate the contribution from the 55 meV peak and the continuum. The total coupling 
constant at T = 290 K from the histogram method is Xhg = 1-85 while Xmfl = 1-9 and Xsf 
= 1.8. In the latter two cases the peak contributions are Xpeak,MFL = 0.88 and Xpeak,SF = 
0.65 respectively. The peak energy is somewhat lower in energy than the low energy dispersion 
kink seen in ARPES experiments [?] . The coupling constants obtained by ARPES are of the 
order A ~ 0.3 - 0.5 [?, I22l ?, ?, [23] as are the coupling constants derived from LDA [?, ?]. 
The remaining coupling constant Xcont ~ 1-2 for the SF contribution also corresponds well with 
earlier estimates [?]. When the temperature is decreased the total coupling constant increases: 
X200K = 2.0 and Aioo_ft: = 2.3. Figure [3}: indicates that this increase arises not simply from an 
increase in coupling to the mode but rather from an increase in intensity in the energy range 
below 50 meV. Since we cannot exclude pseudogap effects to play a role here we refrain from 
separating the contributions due to the mode and the background. 

4. Conclusion 

We have presented a detailed analysis of optical spectra of the high temperature superconductor 
HgBa2Cu04+5 in terms of strong coupling theory. A comparison of the marginal Fermi liquid 
model and the spin fluctuation model shows that the latter better describes the optical spectra 
of Hg-1201. However, the best description of the optical data is obtained when a low energy 
peak is added to either of the two models. A histogram model of n(tj) is dominated by a 
mode at 55 meV and this is lacking in the previous models. The addition of an extra mode to 
the MFL and SF models greatly improves the agreement with experiment for both models and 
allows us to separate the n(a;) spectrum into two components. The first is a nearly temperature 
independent peak centered at 55 meV with a coupling strength at room temperature Xmode ~ 




(5) 



0.7. The energy and coupling strength of this mode strongly suggest an interpretation in terms 
of electrons coupling to a vibrational mode. The second feature is a continuum with a maximum 
around 200 meV with a coupling strength of about Xcont ~ 1-2. 
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Appendix A. 

The calculation of the optical conductivity based on n(u;) rests on several assumptions. The 
most important are (i) the system is a Fermi liquid with a constant density of states, (ii) the 
electrons are coupled to bosonic modes that have little or no dispersion in k-space and (iii) one 
can ignore vertex corrections. For conventional superconductors these assumptions have been 
shown to be reasonable but they are not apriori correct when strong correlations effects are 
important. Calculations have shown that these assumptions still apply when the bosons are 
spin fluctuations with a spectrum that extends to relatively high energy in the strong coupling 
regime [20j. With the above approximations the complex optical conductivity can be 

expressed in terms of a self energy in the following way |24| [25] . 



+00 



npiu) + X, T) — nF{x, T)dx 
iA-KUJ J oj — S(x + w, T) + S*(x, T) + iVimp ' 



(A.l) 



where ujp is the plasma frequency, Timp is an impurity scattering rate, npix) = {ex p{l3x ) + 1) ^ 



is the Fermi-Dirac distribution function and (3 = {k\)T) ^ . The self energy in Eq. A.l 
by, 

OJ — e + uj' + i5 OJ — e — oj' — i5 



E (cj, T)= de doj'li{oj') 
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(A.2) 



where nB{x) = {exp{(5x) — 1) ""^ is the Bose-Einstein distribution function. The integral over e 
can be performed analytically, leaving. 



S(w,r) = / doj'Ii{J)L 
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(A.3) 
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(A.4) 



The in this last expression are digamma functions. Equations (A.l A.4) together with the 
standard Fresnel equations allow us to calculate any optical property given n(a;) . From these 
equations one can also immediately see that it will be very difficult to extract n(u;) directly from 
the experimental data: this would require at least a second derivative of the optical data . 
Given this difficulty we use the following approach. We choose an analytic form for n(u;) with 



several adjustable parameters and numerically perform the integrations in equations (A.l) and 



(A.3). The parameters of n(a;) are optimized using a standard Levenberg-Marquardt algorithm. 
This algorithm revolves around minimizing the function. 
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(A.5) 



where the yi is the experimental datapoint value corresponding to the point Xi and 
f{xi,pi,...,pm) is the calculated value in the point Xi from parameters pi...pm that are to be 
optimized and the are the error bars on the value yi. Since the optimization process requires 
knowledge of error bars on the quantity to be fitted, we use the reflectivity spectra measured in 
the infrared region of the spectrum for which we have accurate estimates of the error bars in- 
volved. Note that the calculation of the conductivity requires a choice for the plasma frequency, 
tOp. We add iVp as an independent fit parameter, but check that it is consistent with values 
obtained from a spectral weight analysis. The method described above applies to the intraband 
conductivity only. However, our spectra contain contributions due to interband transitions as 
well [9]. These we model using standard Lorentz oscillators added to the conductivity calculated 
using equation (A.l). A second point to be noted is that our analysis can only be carried out 
in the normal state. The analysis of spectra in the superconducting state is more challenging 
since it requires solving the full Eliashberg equations with high precision. 
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